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Abstract 

A recently developed formulation for a direct treatment of the equations for two- and three-nucleon 
bound states as set of coupled equations of scalar functions depending only on vector momenta is extended 
to three-nucleon scattering. Starting from the spin-momentum dependence occurring as scalar products in 
two- and three-nucleon forces together with other scalar functions, we present the Faddecv multiple scattering 
series in which order by order the spin-degrees can be treated analytically leading to 3D integrations over 
scalar functions depending on momentum vectors only. Such formulation is especially important in view of 
awaiting extension of 3N Faddeev calculations to projectile energies above the pion production threshold 
and applications of chiral perturbation theory 3N forces, which are to be most efficiently treated directly in 
such three-dimensional formulation without having to expand these forces into a partial wave basis. 



PACS numbers: 21.45.-v, 21.30.-x, 25.10.+S 



I. INTRODUCTION 



The three-nucleon (3N) system is the first nontrivial case to learn about the action of nucleon- 
nucleon (NN) and three-nucleon (3N) forces in bound states and scattering observables. Below 
about 200 MeV laboratory projectile energy 3N scattering can be well treated in a momentum 
space representation based on the Faddeev equations in a partial wave 

elastic nucleon-deuteron (Nd) scattering as well as breakup processes. A different approach using 
a coordinate space representation and based on hyper spherical harmonics is equally precise for 
elastic Nd scattering [2j. 

Over many years a rich set of Nd data has been accumulated , which is not only a 

very valuable source of information about the spin- and momentum-dependence of nuclear forces 
but also about the reaction mechanism of multiple rescattering processes. The nuclear forces un- 
der consideration today are on the one hand the so-called (semi)-phenomenological high precision 

n □ a 

forces [a, 13, 13], describing the NN data up to the pion production threshold perfectly well. The 
three-nucleon forces on a semi-phenomenological level are much less developed |9j and are not 
constructed in a consistent manner with respect to corresponding NN forces. Nevertheless those 
NN and 3N forces describe 3N data (elastic as well as breakup cross sections and numerous spin 
observables) often spectacularly well. However, there are exceptions where one finds serious dis- 
crepancies between that theoretical prediction and the data, especially for some spin observables 
at the higher region of that energy range 3, J, 5|. 

In recent years effective chiral perturbation theory links nuclear forces to the symmetry of QCD 
and diagrammatically builds up the nuclear forces in a systematic expansion [10]. In this fashion 
NN, 3N and even 4N forces are at present consistently generated up to next-to-next-to-next lea ding 
order (N 3 LO). Applications of these forces in the few-nucleon sector [111 ] and for light nuclei [12j 
are quite successful. However, the effective theory is limited in the energy regime it can be applied 
to, which is related to the smallness of the parameter underlying that expansion. A typical upper 
limit for the applicability of those chiral forces is 100-150 MeV nucleon laboratory projectile energy. 

Considering the nucleon projectile energies above the pion production threshold and even going 
into the GeV region, no similar systematic approach to nuclear forces is yet available. However, 
this energy region imposes challenging questions not only about the nuclear force, but also about 
the underlying reaction mechanism. One should expect an increased importance of 3N forces , the 
influence of baryon resonances, meson production, to name a few. Another challenging question 
is an investigation of the limit, where hadronic degrees of freedom apply and where subnuclear 
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degrees of freedom must be explicitly considered. In order to enter this energy region in the 3N 
sector various challenges have to be overcome. On the technical side the standard partial-wave 
decomposition (PWD) has to be given up due to the strongly increasing number of partial wave 
states that need to be summed for a converged result. To face this challenge the direct use of 
momentum vectors, i.e. a 3-dimensional (3D) formulation of the problem turned out to be a 
promising path. This has been documented in various studies for three-boson scattering carried 
out in a Faddeev scheme in momentum space [13]. In addition, the high energy region requires that 
Galilean invariance has to be replaced by Poincare invariance. For the case of three-boson scattering 
this has successfully been achieved fl^l. fis|. where the Poincare invariant Faddeev equations have 
been solved for projectile energies in the GeV regime. As already said the enormous challenge will 
be to develop the underlying dynamics. 

It turns out that Poincare invariant formulation of the 3N system is already important at quite 
low laboratory energy. Relativistic effects are discernable in some regions of the breakup phase- 
space starting at the energy of the incoming nucleon about 65 MeV [lq |. For Nd elastic scattering 
vector analyzing power they contribute to the famous analyzing power puzzle at energies around 
10 MeV [13]. 



In this paper we shall focus on the question how to incorporate spin and isospin degrees of 
freedom into the bosonic three-body calculations so that the successful approach of Ref. [r| is 
applicable to nucleons. The aim is to reduce the formulation with spin/isospin degrees of freedom 
to scalar, spin independent functions of vector momenta in the same spirit as already presented 
for the 2N and 3N bound states in Ref. [18]. The idea is to use the original structure of the NN 
forces consisting of scalar operators in spin- and momentum-space and scalar functions which only 
depend on momenta. This carries over to the NN t-operator which is a central building block in 
the Faddeev scheme. In addition, 3N forces appear naturally in this formulation. 

One form of the Faddeev equations for 3N scattering is based on the multiple scattering series 
for the breakup process which can be summed into a Faddeev integral equation for a transition 
operator T|$) in our standard notation [19|]. Here |<E>) is the initial product state of a deuteron 
and a momentum eigenstate of the projectile nucleon. Analogous to the NN t-operator there will 
be an operator form for three-nucleon transition operator T|<I>). However the number of scalar 
spin-momentum operators will be enormously high, which makes it not advisable to rewrite the 
Faddeev equation for T|$) into a coupled set of equations for the accompanying scalar momentum 



dependent functions. The experience with high energy three-boson scattering 



131 ] suggests that it 



is promising to generate instead the multiple scattering series, which will automatically generate 
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the operator expansion order by order. We want to use this insight as starting point for our study. 

In Section [TT] we introduce the necessary formal ingredients. The lowest order in the multiple 
scattering series is worked out in Section IIII1 The next term, which is second order in the NN 
t-operator, is constructed in Section IIV1 In this order, the free 3N propagator appears for the first 
time, leading to the notorious moving logarithmic singularities. This can be avoided totally as 



shown in 
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2l| . We apply this new method in which the logarithmic singularities are replaced 



by single poles, which then can be handled in a similar fashion as the poles in the 2N Lippmann- 
Schwinger equation. Having the second order under control it is obvious to go to the next order. 
However, we will not give this obvious continuation explicitly. The inclusion of 3NF's, which is 
most interesting from a physical point of view will be skipped in this first formal attempt. However, 
we do not expect principal difficulty according to the experience for the 3N bound state, where 
the inclusion of 3NF's has been worked out explicitly [18]. The calculation of observables based on 
the 3-dimensional form is described in Section [V] Various Appendices provide further information. 
Finally we conclude in Section IVll 



II. THE FORMAL INGREDIENTS 

Our standard form of the Faddeev equation for 3N scattering is given by 

T\4>) = tP\<t>) + (1 + tG Q )vi l) (1 + P)\4>) + tPG T\4>) 

+ (l + tG o )y 4 (1) (l + P)T|0) , (1) 

where t is the t-operator for the NN pair 23, Go the free 3N propagator, P the sum of a cyclical 
and an anticyclical permutation and vjj^ the part of the 3NF which is symmetrical under exchange 
of nucleons 2 and 3. Here we arbitrarily choose nucleon 1 as being the spectator. 

Knowing T\(j>) amplitudes for Nd elastic and breakup scattering are given by the matrix ele- 
ments [l|, H ] 

($'|[/|$) = (^PG^ 1 + PT|$) , (2) 

<$o|tfo|$> = <$o|(l + P)T|$> . (3) 

The state describing three free nucleons is given by | <£())• 

The iteration of Eq. (pQ) generates the Faddeev multiple scattering series. When neglecting 
3NFs, the iteration leads to 

T|$) = tP\<S>) + tPG tP\<S>) + ■ ■ ■ (4) 



As we learned from 3-boson scattering driven by a 2-body force of Malfliet-Tjion type [l3|] which 
incorporates typical properties of the NN force, namely an intermediate range attraction and a short 
range repulsion, the convergence of the series from Eq. improves with increasing energy [^J. 
First, we introduce the three possible 3N isospin states |7 a ) = |(i a ^)T a Mr): 

bo) = \(q\)\m t ) , 
| 71 ) = \(i\)\m t ) , 

| 72 ) = |(1~)~M T ) , (5) 
in which the 2N isospin t is coupled with the isospin \ of the third particle to the total isospin 



T = \ or f . As is well known 



241 ] in isospin space the 2N t-operator has the form 



t = Yl l7a)WTb| • ( 6 ) 
ab 

We assume conservation of t a but allow for charge independence and charge symmetry breaking 
which leads to the coupling of T = \ and | states: 

tab = $t a t b tt a T a T b ■ (7) 

The linear combination of np t-operators in t = and t = 1 states with pp (nn) t-operators 
for proton-deuteron (pd) (neutron-deuteron (nd)) scattering is given in [25f] . Furthermore, the 
permutation operator P in the 3N isospin space reads [3] 

<7a|P|76> = S TaTb F taTaTb (P^PiP + (-Y^P^P^) , (8) 



where Ft a T a T b is essentially a 6j-symbol 24| and P?™ are transpositions of the nucleons ij acting 



only in spin and momentum spaces. The 3N momentum space is spanned by the standard Jacobi 
momenta p and q [z]]. Combining spin/isospin space with the momentum space leads to the 
permutation operator 

5{p- rfffi'^P&Ph + (-) ta+tb S(p + *(qq '))8(p ' + n'(qq'))P? 3 Pi 3 (9) 



with 



Next we use the operator form expansion of the off-shell NN t-operator [ieI in the 3N spin 



and momentum spaces 

(ffltt a T a T b \p'q') = ^t^ TaTb (p,p\E q ) Wj (a(2),a(3),p,p')S(q-q') , 



(11) 



where E q = E — ^q 2 and E the total cm. energy. Due to parity and time reversal invariance 



n 



exactly 6 terms of scalar spin-momentum operators Wj are possible [27(. They are 



wi(a(i),a(2),p,p 
w 2 (a(l),a(2),p,p 
w 3 (a(l),a(2),p,p 
Wi(a(l),a(2),p,p 
w$(B{l),d(2),p,p 
w e {a(l),a(2),p,p 



<?(!) • <?(2) 

(<?(l) + a(2))-(px#0 
a(l)-(pxp')a(2)-(pxp') 
a(l).(p + p')a(2)-(p + p') 
a(l)-(p-p')a(2)-(p-p'). 

(3) 



(12) 



Using this operator representation of the NN force, the functions t t T Tb (p,p',E q ) of Eq. (jlip are 



scalar functions and depend only on three variables, the magnitudes of the vectors p and p and 
the angle between them, p ■ p ' . 

Finally, we use the operator form of the initial state as given in Ref. 28] 



(p§1(7o|$) = (pq\(p) 

= (01 (p) +Mp)(&( 2 ) -P - -^P 2 )\^m d )\m 10 )5(q - q ) 

2 

= ^(/)k(p)Ok(v(2),a(3),p)\lm d )\mi )5(q - q ) . 
k=l 



(13) 



Here (pi (p) and <f>2 (p) are proportional to the standard s- and d-wave components of the deuteron 
wave function. The state is the pure spin 1 two-nucleon state with spin magnetic quantum 

number m^, and mio is the initial state spin magnetic quantum number of the projectile nucleon. 
The vector % is the initial relative momentum in the Nd system. All these are the necessary 
ingredients to work out the terms in the multiple scattering series of Eq. (J4j) . 

III. THE FIRST ORDER IN t 



We consider the first order term in Eq. (|4|) and use Eqs. 



TT1) and (I13H to obtain 



(pq\(la\tPm = F tatoTo (pq\t taTaTo (P??Pir + (" 



\t a +to r>sm r>sm 



(14) 
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Here we used the isospin property of the initial state: to = 0, To = 1/2. Then 
(pq\(la\tPm 

= F tM f d 3 p'd 3 q'^t^ TaTo (p,p\E q ) Wj (a(2),a(3),p,p 
J 

J d 3 p"d 3 q"(S(p' - n(q'q"))S(p" - 7? '(q " 'q "'))P[ 2 P^ 

2 

^4(/)O fc (?(2),?(3),p ")|lm d )| mi0 )%" - 90) 

k=l 

= FtatoTo E ttaT a T (P> ^(99o) , (<?(2) , (3) , p, 7? (g?o))Ti 5 2 P| 3 



E^(|vT / (^o)|)Ofc(d ! (2),<T(3),7r / (^o))|l"id)|"iio} 
fc=i 

+ (-) ta+ *° E *2t o T (P. #>;(*(2), <?(3), jf, -ff(TO ))A S 3^23 



2 

^^(^'(^DOfc^),^),^^))!^)^^)] . (15) 
fc=i 

Note that we used the fact that the two operators from Eq. (|13p depend quadratically on the 
momenta. Then we define 

Wj (e?(2), a{3), p, n)P{ 2 P^O k (a(2) , a(3) , vf') 
= ^(<t(2), a (3), p, 7f)O fc (a(3), a(l), vfOP^Plg 

= a ifc (a(l),a(2),o ! (3),p;7r,7r / )P 1 S 2P2 S 3 > (16) 

and 

Wj(a(2), a(3),p, -n)Pl 3 Pi 3 O k (a(2),a(3),ir') 
= Wj (<r(2) , a (3) , -vT)O fc (*(1), *(2) , J)P( 3 P^ 

= b jk (*(l),a(2),ff(3),p,i?,7r')P{ 3 P£ 3 . (17) 

The scalar expressions a,j k and have to be worked out such that each a(i) occurs only once. 
The explicit expressions are given in Appendix [XI 

Inserting then Eqs. (fT6|) and (fTT|) into Eq. (fT5|) yields 

(p5K 7o |tP|$) 
jk 

a jk (a(l),a(2),a(3),p, vr, ^0^12^23 1 lm d> l m io 
+ E &7b(P> ^)^(I^D^oTo(-) ta+ *° 



b jk (a(l),a(2),a(3),p,K,7r')P( 3 Pi 3 \lm d )\m 10 ) . (18) 

We see the expected structure, a sum over the product of scalar operators (aj k ,bj k ) multiplied by 
scalar functions. Note that 

bjk = b jk (a(l),a(2),a(3),p,n(qq ),n\qq )) . (19) 

It remains to display the singularity structure for the calculation of the physical amplitudes and 
the treatment of the second order. 

The NN t-matrix in the pd isospin space has the general structure 

l taT a T h = S t a 0t o li + S ta lt 1Ta T b , (20) 

where 

t ii = , (21) 
tlT a T b = OLT a T b tn p + ^TaT^pp 1 , (22) 
and aT a T b , Pr a T b are numerical values related to Clebsch Gordon coefficients. 



The np t-matrix t^L for t = has a deuteron pole 



t° n ° p (p,p',E q ) 
E g + ie- E, 



Therefore, Eq. (fTBj) can be decomposed further, 



(pq\(la\tP\§) 



E q + ie - E d 



a jk (a(l),a(2),a(3),p, vf, jf )Pf 2 P^\lmd)\m 10 ) 



yi? E q + ie - E d 
b jk (a(l),a(2),a(3),p,n,n')Pl 3 P^\lmd)\m w ) , (24) 



which retains the structure of scalar spin-momentum dependent operators and momentum depen- 
dent scalar functions. 



IV. THE SECOND ORDER IN t 

Next we consider the second order term in Eq. @ and perform the isospin projections according 
to Refs. ® and © 

{ la \tPG Q tP\<5>) = t ab ( lb \P\ lc )G {j c \tP\$) 
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Ft a t c T b tt aTa T b (P?TPir + {-) ta+t cP^P^)G Q 5 TbTc ( lc \tP\$>) . (25) 



Then we insert Eqs. (fTT|) and ([9]) 



(j><fl(j a \tPG tP\$) 

= Ft a t c T b J d i p'Y J iw b ^v\E q )w j (a(2),a(3),p,p') 

j 

J d 3 p 2 d 3 q 2 (p'q](P{™P£™ + (-)* a+tc P 1 T^T)|P2g 2 )G (p292) 

*zyr c <Rj&K7c|*-P|$) 



n(qq2))5(p 2 - ^ {qq 2 ))P( 2 Pk 

+ (_)*.+*c J(f / + 7f(gT2))^P2 + l?m))P^Piz) 
G (p2q2)ST b T c {P2Q2\{lc\tP\^} 

= F t a t c T b Y^ / t( t a T a T b (P^(^2), E q ) Wj {a{2), a{3),p,jr(qq 2 )) 

Pt 2 PLG^\mq2)6T b T^\m 2 )q2\{ic\tPm 

+ (-)* a+ * c 4fT a T 6 (P> 3>j(*(2), <?(3),£ -ff(gg 2 )) 

A S 3P2 S 3Go(^(^2)g2)^T c (-^(gg2)92|(7c|iP|*)] • (26) 
Finally we use the expression from Eq. (|18|) for the first order term 
(pq]( 7a \tPG tP\<S>) 



= PtatcnY^ J 4kT 6 (^^(^2),^g)^i('?(2),0 ; (3),p,^(^'2)) 
3 

Pf 2 Pi 3 Go(Tr'(qq 2 )q 2 )5 Tb T c 

^2(%T b T if (992) , A02,qo), Eq 2 )Mf (<72 , ft) ) D^ietoTo 
Ik 

ai k (a(l),a(2) 1 a(3),7f'(qq 2 ),7f(q 2 , qo),7r'(q 2 , q ))P{ 2 P 23 \lm d ) \m 10 ) 

+ * £^ bTo ( 9^2 ) , - ^(92 , 0*0 ) , ) 0fc ( I ^ (g*2 , ^0 ) I ) ^t c t „ To ( - ) * C + *° 

blki&Q-), <?( 2 ), 0*(3), ^ / (OT2),^(g2?o) ) ^'(Q2go))A S 3- P 2 S 3l lm <i)l m 'lo)) 

+ (-) ta+tc &T 5 (p> -Ttm),E q ) Wj {a(2), a®,?, -nm) 

Pf 3 Pi 3 G (n\qq 2 )q 2 )5 Tb T c 

(* * ' Wo ( ~ (992 ) , (92 , ft) ) , £g 2 ) 4>k ( 1 7?' ( ft> , 9*0 ) I ) F t c i T 

Zfc 

oijfc^l), a(2),a(3),-7f / (^* 2 ),vr(«2, ft)), ^'(^2, gb))-Pi2- F 23l lm <i}l m io) 
+ ^ ) T 6 To(- ff '(992), -7r(g2, ft), £ g2 )0fc(|7r'(g2, ft))|)-Ft c i To(-) tc+ *° 
^fe(^(l),^(2),CT(3),-w'(gg 2 ),^(g2^),7f'(g2ft))) 
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Pf 3 P^\lm d )\m w )) 



(27) 



Now we collect the spin parts and define 



Cji k (a(l), 3(2), 3(3), p, Tr(qq 2 ), n'(qq 2 ), n(q 2 q ), 7r / (?2?o))^i2- f 23 
= Wj(a(2), 3(3), p, n(qq 2 ))Px 2 P2 3 ai k (3(l), a (2), a (3), Tr'(qq 2 ), T?'(q 2 qo)) 
= Wj (3(2),3(3),p, i((qq 2 ))a lk (3(2), 3(3), 3(l),7f\qq 2 ),n(qlql),7f'(q 2 q ))P^Pi 3 , (28) 

D jik (3(l), a (2), 3(3), p, n(qq 2 ), n (qq 2 ),jf(q 2 q ), n (&q ))Pf 2 P£ 3 
= Wj(a(2),a(3),p, n(qq 2 ))P{ 2 P 23 bi k (3(l), 3(2), a (3), 7f'(g? 2 ), n{q 2 qo), 7r'(g2go)) 
= Wj (3(2),3(3),p, n(qq 2 ))b lk (3(2),3(3), 3(l),Tr'(qq 2 ), n(q 2 q ), n'(q 2 q ))Pf 2 P^ , (29) 

Eji k [a(l), 3(2), 3(3), p, n(qq 2 ), n'(qq 2 ), n(q 2 q ), ^'(g^go))^^ 
= Wj(a{2),a(3),p, -^(^2))A S 3 P 23°;fc(^( 1 )^( 2 )^( 3 ) ; -^'(^2), ^fego), ^'tego)) 
= Wj (a(2),a(3),p,-K(qq 2 ))a lk (aW , (30) 

F^ct(1),ct(2),ct(3),p,tt(^ 
= Wj (3(2),3(3),p, -n(qq 2 ))P? 3 Pi 3 b lk (3(l),3(2),3(3), -n'(qq 2 ), n(q 2 q ), 7r'(q 2 q )) 
= Wj (3(2),3(3),p, -n(qq 2 ))b lk (3(3),3(l), 3(2), -n'(qq 2 ), n(q 2 q ), n' (q 2 q ))P[ 3 Pi 3 . (31) 



Cji k (3(l),3(2),a(3),p,n(qq 2 ),n'(qq 2 ),n(q 2 q ^ 
+ (^*2 ) , (g*2 , ^0 ) , ^ 2 ) ^fc ( I ^ (g*2 , ^0 ) I ) toTo ( - ) ' C +t ° 

D jlk (3(l), 3(2), 3(3), p, n(q%),^\q^ 2 ),n(q 2 ql),n\q 2 q ))) 
+ (-) ta+ic t!f TQT6 (p,-vf(^ 2 ),^)Go(7r'(^2)g2)^T c 

Yl ^tcT b T (-^(992) , AQ2 ,%), E q2 )<t> k ( |7f'(g 2 , 9o) |)*Wzb 



Inserting these expressions into Eq. (|27p yields 




+ t. 



Eji k (3(l), 3(2), 3(3), p, n(qq 2 ),n\qq 2 ),n(q 2 q ),Tr'(q 2 qo)) 

4?T fe T (-^ (99*2) , - 7T (g 2 , ft)) , E q2 ) 4> k ( j 7?' g ) I J-PfetdTo (-) 
Eji k (3(l), 3(2), 3(3), p, n(qq 2 ),nXqq 2 ),n(q 2 q ),n\q 2 q )) 
Pf 2 Pi 3 )]\lm d )\m 10 ) . 



(32) 
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As will be shown for some examples in Appendix B the coefficients from Eq. (|28p together 
with the other coefficients from Eqs. (|29|) through (|31|) will depend on the following types of scalars. 
As an aside, one more term, a{i\) ■ q 2 (&(i 2 ) x " g*2, is expected to occur in Cji k beyond the 

ones given in Appendix C: 

Cji k (a(l), a(2), a(3),p, Tr(qq 2 ), n'(qq 2 ), n(q 2 q ), t?'(<Mo)) 
= C$ (<?(1), <x(2), (7(3), p, q, q ) + C®(p, q, q , q 2 ) 
+ ^Cj^i^^^^^P^^o^ih) -qi 

h 

+ Yj C i?fc;iii 2 ^^ 3 )'^' & 9o)^(il) ■ q2#(h) ■ <?2 

U«2 

+ Cjifc^' <fo)<?(n) • <72<?(*2) • gVfe) • <7*2 

+ S C jilht 2 (°( i3 ^Pi^ 9bM*i) x <?fe) • 42 

+ E C jffc;n i2 J£?>9o)^^ (33) 

where 11,12,13 = 1,2,3 and cyclical permutations thereof. 

It is sufficient to consider the first integral in Eq. ([27|) going along with the coefficients Cjik- 
We insert Eq. ([33]) : 

(Ml( 7a |tPG tP|$) (1) 

/ d3 Q2t ( t J a ) TaTb (p, T?(qq 2 ), E q )G (jr' (qq 2 )q 2 ) 
t( t% b T (n\q&),n(q 2 ,ql),E g2 )<p k (\n\q 2 ,qo)\) 
+ J cPq 2 t { £ TaTb {p, K(qq 2 ),E q )G (n(qq 2 )q 2 ) 

^TsTo^'^),^®, go), E g2 )(f> k (\Tf'(q 2 , gb)l)cj£(p, g, g , g 2 ) 

+ E^ijK^),^),^^^) / ^ 3 g2^T a T i ,(P^(^2),^ (? )Go(^ / (9?2)g2) 
1 

, 



t ilkro(^ / (^2),^(g2,go),Pg2)^A : (|v? / (g2,go)l)o ! (^i) ■ h 



+ E^iffeji^^fe),^,^^) y d 3 ^2tSf TaT6 (p,^(^2),P g )G (^ / (^2)g2) 

U«2 

t t?T 6 To(^ / (^2),^(g2,go),Pg2)^(|7?'(g2,go)l)<7(n) • g^fe) • <7*2 



+ C^l(p,q,q ) J d 3 g2^kT b (^^(gg2),Pg)Go(7f / (gg 2 )g 2 ) 

tunTo^'i^),^^,^), E q2 )(f) k (\Tr'(q 2 , q )\)ff(ii) ■ q 2 B{i 2 ) ■ q 2 a(i 3 ) ■ q 2 

+ ZJ c ]fci 42 ( <5; ( i 3),P,g"',g*o) y d 3 g24iT t (^^g2),^)Go(^(gg2)?2) 



217^2 

H% b T (K\qq 2 ),T?(q 2 , q ), E q2 )(j) k (\7r'(q 2 ,q )\)a(ii) x a(i 2 ) • g 2 
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+ E C jt,hi2iMZ<l°) f d\ 2 t^ TaTb {p^(qq 2 ),E q )G (7f'(qq 2 )q 2 ) 

4?T 6 T (^V<?2),7r(<72, 9o), E q2 )<p k (\Tr'(q 2 , q )\)ff(ii) ■ q 2 {B{i 2 ) x d*(i 3 )) • q 2 
P{ 3 P£ 3 \lm d )\m 10 ) ■ (34) 

The integrals in Eq. (|34p are of the following form 

*l = J d 3 q 2 Y , (35) 

X 2 = J d 3 q 2 Y a(i) ■ q 2 , (36) 

X 3 = J d 3 q 2 Y a(i) ■ q 2 a(j) ■ q 2 , (37) 

X 4 = J d 3 q 2 Y a(i) ■ q 2 a(j) ■ q 2 a(k) ■ q 2 , (38) 

X 5 = J d 3 q 2 Y a{t) x a(j) ■ q 2 , (39) 

X & = j d 3 q 2 Y a(h) ■ q 2 {B{i 2 ) x a(i 3 )) ■ q 2 , (40) 

where Y represents functions with different dependencies on scalar quantities. For instance in the 
first integral in Eq. (|34p on has 

H% b T (rf\qq2),n(q 2 ,qo),E q2 ) c/)k(\n'(q 2 , q )\) . (41) 

As detailed in the Appendix C its angular dependence is 

Y = Y(p ■ q,p -q 2 ,q- q 2 ,q- q ,q 2 ■ %) ■ (42) 

The remaining integrals will be treated as exemplified for X 2 . The integral has the form 

X 2 = a s (i) J d 3 q 2 Y q 2s . (43) 

If Y does not depend on q 2 then X 2 = 0. Otherwise, due to the scalar nature of Y, one must have 
the following structure 

J d 3 q 2 Y q 2 ^ s = ap s + (3q s + -yq 0s , (44) 

and therefore 

X 2 = aa(i) ■ p + f33(i) ■ q + j<r(i) ■ % ■ (45) 
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The scalars a,j5, 7 are then determined by multiplying Eq. (|44[) by p s , q s , and qo s , respectively. 
This leads to three equations 

J d 3 q 2 Yp-q 2 = ap 2 + (3p- q + jp- q , (46) 
J d 3 q 2 Yq-q 2 = aq ■ p + fiq 2 + jq ■ q , (47) 
J d 3 q 2 Yq ■ q 2 = aq ■ p + f3q ■ q + 7^ , (48) 

for a, f3, 7. The three integrals have to be determined numerically. 
Correspondingly, despite being more involved, X3 can be determined 

X 3 = a s (i)a t (j) J d 3 q 2 Yq 2s q 2t . (49) 

One has three external momenta in the function Y, namely p, q , and qo. Thus, there are at most 
9 possibilities: 

9 



/ d 3 q 2 Y q 2s q 2t = £ AlQlsQ'u , (50) 
J L=l 

where the QlsQ'h belong to the set 

QlsQ'li = [PsPt,p s qt,p s qot,qsPt, q s qt, q s qot, qosPt, qosqt,qo s qot] • (51) 

By appropriate multiplications from the left one will have nine equations for the coefficients Al, 
generated by 9 corresponding integrals. Thus X% takes the form 

9 

X 3 = ^2 A L <?(i) ■ Q L B{j) ■ Q'l • (52) 

L=l 

Furthermore, 

X 4 = a s (i)a t {j)cF u (k) j d 3 q 2 Y q 2s q 2t q 2u ■ (53) 
Now we have 27 possibilities and consequently 27 equations need to be solved with the result 

27 

X, = Y, B L S{i) ■ Q L a(j) ■ Q' L o(k) ■ Ql , (54) 

L=l 

where QlsQ'uQ'lu are ou * °^ the se ^ 

QlsQ'uQ'lu = \PsPtPu,PsPtqu,PsPtqou, ■■■} ■ (55) 
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The appropriate 27 integrals need to be determined numerically. Finally, X§ can be handled like 
X 2 , and Xq like X3. All the integrals occurring in Eq. (|34|) and the ones following from that are 
of the type 



H = J d 3 q 2 t^ TaTb (p,Tf(qq 2 ),E q ) G (k' {qq 2 )q 2 ) 

4^T (^'(^2),^(i?2,go),^g2)^A : (r7r / fe,go)!)/ , (56) 



where / consist of scalar products of the momenta p, q, %, and q 2 . 

We use the decomposition of the NN t-matrix from Eq. (|20p in the pd isospin space and the 
deuteron pole structure from Eq. (|2"3"|) , Then we obtain 



H = J d\ 2 G Q ^\m 2 )q 2 )(j> k {^\q 2 ,q Q )\)f 



Kotfh (P, n(qq 2 ), Eg) + 5 ta it$ aTb (p, n(qq 2 ),E g ) 
^c0*q1 i (n\qq 2 ), Tf(q 2 , qo),E q2 ) + 5 tc itQ bTb (n\qq 2 ) , T?(q 2 , go), E q2 



2 2 

J3 



d 3 q 2 G {tt' (qq 2 )q 2 ) k ( \n' (q 2 , g ) | ) 

' ilrp' J \p,TT(qq 2 ),E q ) {j) ^ 

St «° E +i € -E d ^iTaTtfa n m2), E q ) 



tn P ' \n'(qq 2 ),n(q 2 ,q ),E q2 ) ^ x +w 
E q2 +ie- E d 



>/.u — h 5t c it^ bTb (n\qq 2 ),jT(q 2 ,q ),E q2 ) \ . (57) 



In numerical implementation one can either follow the standard path [1] dealing with the free 
propagator singularity and leading to moving logarithmic singularities, whose treatment is well 
controlled and documented in [l^ or one applies the new way [^(j, [21I which avoids that com- 
plication totally. We suggest the second path, which for the convenience of the reader is detailed 
again in the Appendix D. 

Summarizing, the second order amplitude of Eq. (I32p using all information given above will 
have the structure 

<pq\(j a \tPG tP\$) = ^5/(p, q,q )O f (a(l), a(2),ff(3),p,q,q )P f \lrn d )\mi ) , (58) 

/ 

where Of are spin-momentum dependent scalar operators and Sf only momentum dependent scalar 
functions, all depending just on the external momenta p,q, and go- Further Pf is either l,P( 2 P 2 z 

This form will be the input for the next order in Eq. @, namely tPGo(tPGotP\&)) where the 
kernel tPGo is exactly treated as above and the first order term tP\$) is replaced by that second 
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order form. Obviously this will go on like that and the total amplitude (pq\ ( 7a |T|<3?) will appear in 
the form of the right hand side of Eq. (|58|) . 

It remains to present the calculation of the observables, which is described for the example of 
the breakup cross section in the next section. 

V. THE ND BREAKUP CROSS SECTION 

The full breakup amplitude according to Eq. ([3|), projected on spin states and using Eq. (JSj) is 
given by 

(mi | (m 2 1 (m 3 1 (p§1 ( 7a | (1 + P)T\ $} 
= (mi|(m 2 |(m 3 |(pg|(7 |T|$) + (m 1 \(m 2 \(m 3 \(pq\('y a \PT\$) 
= (m 1 \(m 2 \(m 3 \(pq\('y a \T\§) 

+ T / F wJTMrni\{rn2\(m 3 m(PuP^ + H U+tb P^T)^^!*) • (59) 

b 

The simplest approach is to apply the permutations to the left, using 

<fflP%P% = (-^+^,-P-^1 = (p (1 V (1) | , (60) 
m*®P% = (-^P- \% = <P (2) ? (2) l • (61) 



Therefore, 



{mi\{m 2 \{m 3 \(pq\(j a \(l + P)T\$) 
= {mi\{m 2 \{m 3 \(pq\('y a \T\$) 

+ E F ^^( m lK m 2|(m3|Pa s 2 P|3(p (2) g (2) K7 6 |r|$) 

b 

+ ^F tatbTa 5TM-) ta+t Hrni\(m 2 \(m 3 \Pl 3 P^(p «g {1) \(j b \T\^) . (62) 



b 



Using now the general form of the right hand side of Eq. (|58p one obtains 

(mi | (m 2 1 (m 3 1 (pq\ (7a 1(1 + P)T|$) 

= (mi [ (m 2 |(m 3 |E (lP> Q,<to)Of{a(l), a (2) , <r(3) , p, g, g )P/ 1 lm d ) |mi ) 
/ 

+ E Ft M 6 T a T b (mi | (m 2 1 (m 3 |Pi s 2 P| 3 

b 

E S f (p P), g ( 2 ), $>)O f (B(l), a(2), a(3),p (2) , g* (2) , ?o)P/ |lm d )|m 10 > 
/ 

+ E *Ut>.31, R ta+ '> (^1 1 <m 2 1 (m 3 |Pf 3 P| 3 
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E Sf(p q qo)O f (a(l), 3(2), a(3),p « q « g )P / |lm d }|m 10 ) 
/ 

As shown in Ref. jlj], the density matrix for the final state is given by 

(Pfh = J2 N ^)kiNl , (64) 

kl 

where 

Nij = (Ai\{pq\Uo\<t>)\qo)\*j) , (65) 
is the expression of Eq. (j63|) . where the set of spin magnetic quantum numbers is given by 

|Aj) = |mi)|m2}|m 3 } , 

|Aj) = |mio)|lm d ) , (66) 
Furthermore, pi is the density matrix for the initial state, 



ft 6 



l -Tr{ Pl )Y,S»(Snii (67) 



with [S v ] being the complete set of 2N spin matrices [1]. For an unpolarized initial state one has 

Pi = j:Tr(pi)l 2X 2 x 13x3 • (68) 
b 

With this, the unpolarized breakup cross section is given up to a phase-space factor as 

= ^EEKAil^W)!^}! 2 

i k 



= 6 5^ E ( m i m 2m 3 |(p?lC/o|'/ ) )l9o)|l^ ( i)|"iio) 

mim2in3 rridmio 

(lm d \(m 10 \(qo\((j)\U^\pq)\rn 1 m2m3) . (69) 
Using now the final form in Eq. f)63|) we obtain 



6 

(mim 2 m 3 |o a (a(l),a(2),d ! (3),p,Q,go)-Pa|lw ( i)|"iio) 



~ ^ /a/fl (mim 2 m3|o Q ,P Q P / 30^|mim2m 3 ) . (70) 



o,/3 mim 2 m 3 
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VI. SUMMARY AND CONCLUSIONS 



We extended the recently developed formalism for a new treatment of two- and three-nucleon 
bound states in three dimensions to the realm of nucleon-deuteron scattering. The aim is to 
formulate the momentum space Faddeev equations in such a fashion that the equations are reduced 
to spin independent scalar functions in the same spirit as already worked out for the 2N and 3N 
bound states. We use the original, most general structure of NN forces, which is built out of scalar 
operators of spin- and momentum-vectors plus scalar functions that only depend on momenta. 
This structure carries over to the NN t-operator which is a central building block in the Faddeev 
scheme. Generating the multiple scattering series for the Faddeev equations we arrive at the 
general structure of the multiple scattering terms in the form of scalar spin-momentum dependent 
operators and momentum dependent scalar functions. Since our formulation only depends on 
operator and momentum vectors, there are no intrinsic limitations on the energy range in which 
it can be applied. This is contrast to the partial wave projected form of the momentum space 
Faddeev equations. 

Since three nucleon forces appear naturally in the form of scalar operators in spin- and 
momentum-vectors and scalar functions depending only on momenta, the present formulation 
can be extended to include also 3N forces in a straightforward manner. This feature is especially 
important in view of the chiral perturbation theory approach, where the multitude of 3N forces at 
N 3 LO appears to be most efficiently treated directly in this three-dimensional formulation without 
having to expand these forces into a partial wave basis. 
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APPENDIX A: THE SCALAR EXPRESSIONS a jk AND b jk 

Since the scalar expressions aj k and bjk from Eqs. (|16|) and (|17|) have to be determined only 
once and their number is not too large, we provide all of them: 

an = 1 (Al) 
ai2 = -V + (a(l)-vr')(a(3)-vr') (A2) 
a 2i = 3(2) • «t(3) (A3) 

a 22 = («T(l)-7r')(<?(2).vr') 

- i(a(l) ■ 7f')(a(2) x «t(3)) • vf ' 

- V 2 (a(2).a(3)) (A4) 

a 3 i = (<?(2) + <7(3))-(pxvf) (A5) 
3 

«32 = jqo ■ (p x q) (<?(i) • 7?') 

- z(7f'-7f)(a(l)-7?0(^(3))-p) 

+ i(p-n')(a(l)-Tf')(a(3)-n) 

+ (a(l)-7f')(a(2)-(pX7f))(a(3)-7fO 

- V 2 (a(2)+a(3)).(pxvf) (A6) 

a 41 = <t(2) • (p x 7f)a(3) • (p x vf) (A7) 
3 

«42 = t(# X g) • ft(<?(l) ' tt'M 2 ) • (p x tt) 

- ■ n)(a(l) ■ 7t')ct(2) • (pxjf)(ff(3) -p) 
+ i(p-n)(a(l) -n')a{2) ■ (p x v?)(<t(3) - tt) 

- V 2 ff(2)-(px7f)ff(3)'(px7f) (A8) 
a 5 i = (<?(2) • (p + t?))(c?(3) • (p + vr)) (A9) 

«52 = ((p*+7f)-7f')(<?(l)-7f')(a(2)-(p + 7f)) 

+ i(a(l) ■ 7f')(a(2) • (p + vf))(a(3) • (px vr')) 

+ <|(^(1) • 7x0(^(2) • (p*+7?))(a(3) • (gx ft)) 

- \iT' 2 (e{2) ■ (p-+ 7f))(a(3) ■ (p + vf)) (A10) 
a 61 = ((7(2) • (p - vf))(a(3) • (p- vr)) (All) 
«62 = ((p'-vf)-vf')( ( ?(l)-vf')(a(2)-(p-vr)) 



+ • vf')(<?(2) • (p - vf))(a(3) • (p x vf')) 

- • vf')(a(2) • (p - v7))(a(3) ■ (g x ft)) 

- V 2 (a(2) • (p- vf))(a(3) • (p - vf)) (A12) 
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612 = -\<*' 2 + (a(l) • 7f')(a(2) • vr') 

621 = <?(2) • <?(3) 

6 22 = (a(l)-7f')( ( ?(3)-7f') 

+ •7r , )(o ! (2) x ^(3)) -tt' 

- V(a(2)-a(3)) 

631 = -(c?(2)+a(3)) -px 7? 
3 

^32 = -^(p x q) ■ qo(a(l) • tt') 

+ i(7T ■7f / )(^(l)-^)(^(2) -P) 

- i(p-7?')(a(l) -if){a(2) ■ tt) 

- Kl).vfO(a(2).vf')(a(3)-(pxvf)) 
+ V(a(2) + ( r(3))-(px^ 

641 = <r(2) • (p x 7t)ct(3) ■ (p X 7?) 
3 

^42 = ^(pxg> go(^(l) • tt'M 3 ) • (p x vr) 

- i(vf •7f / )(<?(l) •^ / )(<?( 2 ) -pW) • (px vr) 
+ i(p ■ 7f')(a(l) • 7f')(o ! ( 2 ) • ^M 3 ) ■ (p x 7?) 

- -vt /2 (t(2) • (px tt)ct(3) • (p x vr) 

651 = (a(2)-(p-7r))(a(3)-(p-7f)) 

b 52 = (P~ vf) • 7f'(a(l) • vr / )(<?(3) • (p - 7?)) 

+ i(a(l) • 7f')(a(2) • (px 7f'))(a(3) • (p- vr)) 

- • 7f')(a(2) • (g x q ))(a(3) ■ (p - tt)) 

- ^(?(2)-(?-7?))(?(3)-(?-7?)) 

6ei = (a(2)-(p + 7f))(a(3)-(p + 7f)) 

b 62 = ((p + vr)-7r / )(<?(l)-7f , )(^(3)-(p + 7f)) 
+ ■H , ){ff{2) • (px 7?'))(a(3) • (p + 7f)) 

+ \i{a{l) ■ 7r')(a(2) • (g x &))(*(3) • (p + tt)) 

- ^7r /2 (a(2)-(p + 7f))(a(3)-(p + 7r)) 
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APPENDIX B: EXAMPLES FOR THE COEFFICIENTS C jlk 

Again the coefficients Cji^, Dji^, Ejy. and denned in Eqs. {(2SD- d3IJ) have to be calculated 
only once. In the following we provide some examples, inserting Eqs. (jlOp for tt and tt' and removing 
the double occurrences of the same 

Cm = 1 (Bl) 

Cll2 = -g(<f2 + T^qo) 2 

+ ((7(2) • q 2 )(a(l) ■ q 2 ) + (<x(2) • q 2 )(a(l) ■ ±q ) 

+ (a(2) ■ • q 2 ) + (<?(2) • ^q )(a(l) ■ ±q ) (B2) 

C121 = <?(3) • <?(1) (B3) 

C 122 = (<r(2) • <f 2 )(a(3) • q 2 ) + (<r(2) ■ q 2 )(a(3) ■ ~q ) 

+ (<?(2) • ~g )(<?(3) • q 2 ) + (<?(2) • ^-o)(«t(3) • |gb) 

- i[(tf(2) • <? 2 )(a(3) x a(l)) • + (tf(2) • g 2 )(a(3) x a(l)) • |qb) 
+ (<?(2) • ig )(<?(3) x <r(l)) • + (<?(2) • ^ )(<?(3) X a(l)) ■ 

- 5(& + 5«) 2 W)-?(l)) (B4) 
Can = e?(2) • <?(3) (B5) 
C 2 i2 = - l -{q 2 + l -q )\a{2). a (3)) 

+ (<?(1) • q 2 )(a(3) ■ q 2 ) + (<t(1) ■ q 2 )(a(3) ■ ^q ) 

+ (a(l) ■ ig )(^(3) • q 2 ) + {&{!) • y o )(a(3) ■ ±q ) 

- ■ q 2 )(a(3) x a{2)) ■ q 2 ) + (<r(l) • q 2 )(a(3) x a{2)) • ~«f ) 

+ (<?(!) • |«&)(^(3) x <?(2)) • + (a(l) • ~gb)(<?(3) x a(2)) ■ -ft)] (B6) 
C221 = • <?(2) - ia(3) • e?(l) x <r(2) (B7) 

C222 = -»(?(3) • (<72 + 5«&)(^(2) x «t(1)) • (<f 2 + iqb) 

- \{& + ■ 9{2)) ~ <*(3) • (*(1) x *(2))] 

+ (92 + ^o) 2 [l + <?(!) • <?(3) - B{1) ■ a{2) - a{2) ■ a{3)] 

+ (<?(!) • (ffl + \qoW{2) ■ (q 2 + ift)) - (ef(l) • + ~gb))(tf(3) ■ (<? 2 + 

+ (<?(2) • (& + ^<fo))(<?(3) • (g 2 + ^0)) (B8) 
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APPENDIX C: THE ANGULAR DEPENDENCE OF Y 



In order to see the angular dependence of a term like that given in (|41|) we use 

n(qq2) = ^<? + <?2 (CI) 

7r'(<Mo) = -<?2 - 2<fo (C2) 

7r(929o) = \p2 + 9o (C3) 

7?'(^2) = -9 - ^92 (C4) 



Therefore 



P-Tr(qq 2 ) = , ^, (C5) 

W(qq2)\ = J<i + \<& + 92 ■ 9*0 (C6) 
*//^.x ^ / — * — * \ I?" 92 + 9 -90 + 3^1 + ^92 -go 

vr (992) • vrfego) = — l-iii-- I ^ ( C7 ) 

|9+ 192II292 + 9o| 



We put q = z and define 





X 


= 92-9 




Xp 


= P - g 




■ 92 


— XpX ~\~ ■ 






= 90-9 


92 


• 90 





(C8) 
(C9) 



x 2 cos((/> p - 02 ) (CIO) 

(en) 



+ Vl ~ x 2 ^/l - x 2 qo cos(fo ~ <P qo ) (C12) 
This settles the angular dependencies of Y. 

APPENDIX D: TREATMENT OF SINGULARITIES 

Regarding the expression (|57|) we face two types of integrals, where Go appears alone or together 
with the deuteron pole: 

Hi = j d 3 q 2 g(p,7r(qq}>),7T\qq 2 ),7i(qlql),^ (Dl) 

H2 = J d 3 q 2 h(p, Tf{qq 2 ), ^'(992), ^{qlqo), n'(q 2 q )) 

Gotf * w (D2) 

bj q2 +ie — tid 
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Here g and h are regular scalar functions. In order to arrive at the new way 21[ one has to go 
back and rewrite both integrals into 

fli = J d 3 q 2 J d 3 p' J d 3 p 2 5(p' -n(qq 2 ))6(p 2 -jr\qq 2 )) 

9(p,fi,P2,K(q 2 qlj),n\q2qo))Go(p2,q2) (D3) 

H 2 = J d 3 q 2 J d 3 p' J d 3 p 2 5(p ' - Tr(qq 2 ))5(p 2 - n\qq 2 ))h(p,p\p2,^(q 2 qo),n\q 2 qo)) 

G (p 2 ,q 2 )-—± (D4) 

E q2 + ie - E d 

Then we change the product of the following two (^-functions 

5{p - n(qq 2 ))S(p 2 - n'(qq 2 )) 



8(p' ~ y^<7 2 + q\ + qq2x)5{p 2 - ^ q 2 + ~g| + Q<& x ) 

-^S(x-x )Q(l-\x \) 
QQ2 



with 



S(P2 ~ \l\<? ~ \ql+p' 2 m\q 2 - \ql+A , (D5) 



n' 2 - i« 2 - a 2 

x = 3 . (D6) 

qq2 



We start with Hi and insert (D5) 

Hi = j d 3 q 2 J dp J dp2 J dp'dp 2 ^-6(x - x )6(l - \x \) 



/ 3 3 3 3 

8(P2 - J^q 2 - -^q\ +p' 2 )Q{^q 2 - ^q 2 + p' 2 )g{p,p \p2^{q2qo)^' imo)) 

— : — x - 2 p1 kp' - - #{<m) ■ m 

E ~\~ it 4rrt^2 m 



Then we carry out the P2 , P2 and p' integrations 

Hi = I J dq 2 J dp'p'dq 2 q 2 8(x-xo)®(l-\xa\) 



3 3 

Q ^4 q2 ~ 4^2 +P' 2 )9(p,p'^(qq2),P2^'(qq2),^(q2qo),^\q2qo)) 
1 



E + K-^ql-i{\q 2 -\ql+p' 2 ) 

2 f dp'p' J dq 2 q 2 @(l - |x |)6(^ 2 - \q\ + p' 2 ) 
g{p,p',q,q 2 )-——- — 2 ri2 ' ( D8 ) 

E + ie - - ±p' 



q 



with 



P2 = \l\{q 2 -q 2 2)+p' 2 , (D9) 
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and 

9(P,P',Q,Q2) 

= J dq 2 5(x - x )g(p,p'n(qq 2 ),p2n\qq2),^(q2qo),^(q2qo)) ■ (D10) 

Remember x = q 2 • q, thus it appears natural to choose q = z and then the x-integration can 
be carried out trivially. 

The two 0-functions restrict the integrations in q 2 and p' into an area whose size depends on 
the magnitudes of the spectator momentum q . It results 

2 f°° 1 



Hi = - dp p' 



qj E + it-^f-y* 

fi+P' 

/ dq 2 q 2 g(p,p ,q,q 2 ) • (Dll) 
J\i-P'\ 

The important point is that the free propagator appears now as a simple pole. 

Next comes the ^-integral, where two singular denominators appear. We rewrite it using again 
(D5) and obtain 



H 2 = J d 3 q 2 J dp' J dp 2 J dp'd P 2^5{x - x )Q{l - \x \) 



3 3 3 3 

S(P2 ~ \/^q 2 - -^q\ +P ,2 )@(^q 2 ~ ^q2+P' 2 )HPiP'>P2,K{&%),K\&q )) 

Sip' ~ TT{qq2))5{p2 ~ n'm)) ■ (D12) 



E + ie- ^q 2 2 -%E q2 +ie- E d 



Then we carry out again the p2 , p2 and p' integrations 

2 



Ho = 



dq 2 J dp'p' dq 2 q25{x - x )6(l - \xq\) 



3 3 

@(^ 2 - -^ql + p' 2 )Hp^p^q^q2) 
l l 

E + ie- it'll - P i Eq2+ie ~ Ed ' 



(D13) 



with 



h(p,p',q,q 2 ) 

J dq 2 5(x - x Q )h(p,p'Tc(qq 2 ),P2^\qq2),^{q2qo),^\q2qQ)) • (D14) 



We rewrite 



E + ie-^-^E + te-^qi-E d 
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-E d -i,4 + Up' 2 + b 2 )' (D15) 

The new denominator function 

6(4feP ' )5 ^iiWTFi (D16) 

cannot become singular inside the integration domain p' — q 2 . Using the expression (D9) we see 
that 

G(q,q 2 ,p') = 1 — 2 = . . ~ r~ 2 • ( D17 ) 
~^d + ^P 2 \ h d\ + ^P 2 



Therefore (D13) turns into 



H 2 = - J dq 2 J dp'p'dq 2 q 2 5(x - x )@(l - \x \ 
3 3 

@(q<? 2 - -^2, +p' 2 )Hp,p',q,q 2 ) 
1 1 



l _ 3 2 P - F + if 3„ 2 — W ]G{1^P') (D18) 

E + ie- ^q 2 - 2_ S + *e - - ^ 

We use the integration over the area in the p' — q 2 plane and obtain 

2 f°° 1 ~ ff+p' 

H 2 = - dp'p' #h(p,p',q,q2) / dq 2 q 2 G(q,q 2 ,p') 

2 f°° 1 , fi +q2 ~ 

/ dq 2 q 2 ——. 3—3 —h{p,p,q,q 2 ) dp p G(q,q 2 ,p) . (D19) 

qJo E + ie--^q 2 2 -E d J\i_ q2 \ 

In the second integral we integrated first over q 2 and then over p'. In both cases we have just a 
simple pole either in p' or in q 2 . 
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